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We discuss the implications of assuming a four-zero Yukawa texture and a general Higgs potential 
for the production of neutral Higgs boson pairs at 77 colliders through the 77 — >■ (fyicf)] {(jji = h, H, 
A) reaction within the context of the two Higgs doublet model type HI. Exact analytical expressions 
for the 77 — > 4>i4>j reaction are presented. The use of a nonlinear iij-gauge, which considerably 
simplifies the loop calculations and renders compact analytical expressions, is stressed. We show 
that these processes are very sensitive to a general structure of the Higgs potential that impact the 
triple and quartic couplings of the scalar sector. We present results for scenarios of the parameters 
of the model that are still consistent with current experimental constraints. It is found that the 
cross sections for the 77 — > (piXpj processes can be up to two orders of magnitude larger than those 
gotten in 2HDM type I and type II. The possibility of a light CP-scalar is also studied. 

PACS numbers: 12.60.Fr, 14.80.Cp 



I. INTRODUCTION 



The discovery of a Higgs boson (or several Higgs bosons) is central to the broad experimental programs of both the 
Large Hadron CoUer (LHC) and the International Linear Collider (ILC) l], either through e^e+ collisions or via the 
secondary 77 and modes 0, IS*] . If there is a Higgs boson, it is almost certain to be found at the LHC and its mass 
measured by the ATLAS [3] and CMS [5|| experiments. Then, in the much cleaner environment of the ILC a more 
complete and precise experimental analysis can be carried out to verify if it corresponds to a Standard Model (SM) 
Higgs or other type of scalar particle. The high level of complementarity between both type of colliders has broadly 
been studied by diverse groups IS]. Whereas LHC has a large mass reach for direct discoveries due its high collision 
energy, the ILC enables precise measurements and therefore detailed studies of direct productions of new particles as 
well as high sensitivity to indirect effects of heavier new particles. It has been found that some processes occurring 
via 76 or 77 collisions are complementary to their analogous reactions at e+e^ collisions as the former are more 
appropriate to study the bosonic sector of the SM. In particular, 77 collisions have been recognized as an invaluable 
tool to probe the structure of electroweak interactions at high energies, both in the gauge and the Higgs sectors [7|. 
In the gauge sector, the reaction 77 — WW has been widely studied, mainly for testing any physics beyond the SM 
[1, As for the Higgs sector, once a Higgs boson is discovered, the one- loop process 77 — )■ — AT might play 
an important role in determining some properties of this elusive particle, such as mass, total width, CP properties, 
and couplings to other light particles in a model-independent fashion [l^]. Recently, a series of papers have been 
published in photon-photon collisions on the more traditional types I and II of 2HDM, which spelt out the genuine 
phenomenological features that differentiate them from the MSSM [ll| . In addition, neutral particle pair production 
at a 77 collider can be highly sensitive to new physics effects as processes of this kind are naturally suppressed because 
they first arise at the one-loop level, thereby providing a detailed test for the structure of extended Higgs sectors. 

It is expected that the LHC will allow us to test the mechanism of Electro- Weak Symmetry Breaking (EWSB), which 
represents a unique probe of a weakly-interacting theory, as is the case of the Minimal Supersymmetric Standard Model 
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(MSSM) [12^ and general 2HDMs of Type I, II, X and Y (2HDM-I, 2HDM-II, 2HDM-X and 2HDM-Y) 

or whether strongly- interacting scenarios are instead realized, like in the old Technicolor models or the ones discussed 

more recently [l^. 

The 2HDM-II has been quite attractive to date, in part because it coincides with the Higgs sector of the MSSM, 
wherein each Higgs doublet couples to the u- or c?-type fermions separately. However, this is only valid at tree-level 
[20L HH- When radiative effects are included, it turns out that the MSSM Higgs sector corresponds to the most 
general version of the 2HDM, namely the 2HDM-III, whereby both Higgs fields couple to both quarks and leptons. 
Thus, we can consider the 2HDM-III as a generic description of physics at a higher scale (of order TeV or maybe 
even higher), whose low energy imprints are reflected in the Yukawa coupling structure. The general 2HDM has a 
potential problem with flavor changing neutral currents (FCNC) mediated by the Higgs boson, which arise when each 
type of quark {u and d) is allowed to couple to both Higgs doublets, and FCNC could be induced at large rates that 
may jeopardize the model. The possible solution to this problem of the 2HDM involves an assumption about the 
Yukawa structure of the model. Then, in order to keep the FCNC problem under control, we can choose one of the 
following mechanisms: (1) Discrete symmetries. This choice is based on the Glashow-Weinberg's theorem concerning 
FCNC's in models with several Higgs doublets [22*1. (2) Radiative suppression. When a given set of Yukawa matrices 
are present at tree-level, but the other ones arise only as a radiative effect. This occurs for instance in the MSSM, 
where the type-II 2HDM structure is not protected by some symmetry, and it is transformed into a type-Ill 2HDM 
through loops-effects of sfermions and gauginos ^2^ [2J, . (3) Flavor symmetries. Suppression of FCNC effects can 
also be achieved when a certain form of the Yukawa matrices that reproduce the observed fermion masses and mixing 
angles is implemented in the model, which is then named as THDM-III. This could be done either by implementing 
the Frogart-Nielsen mechanism to generate the fermion mass hierarchies [2^ , or by studying a certain ansatz for the 
fermion mass matrices ^3- It should be noted, that when this scheme is implemented, the scalar Higgs potential 
must be expressed in its general form because it is not necessary to impose a discrete symmetry. 

In this paper, we will work in the context of the mechanism (3) by implementing the so-called four-texture ansatz 
for Yukawa matrices. This ansazt, is illustrated in Refs. [2^, |2^ [s^], in which a detailed study of the 2HDM- 
III Yukawa Lagrangian was presented under the assumption of a specific texture pattern [sH, which generalizes 
the original model of Ref. 32]. Phenomenological implications of the neutral Higgs sector of the model, including 
Lepton Flavour Violation (LFV) and/or Flavour Changing Neutral Currents (FCNCs) effects, have been studied 
previously (25l . Issj . The extension of such an approach to investigate the charged Higgs boson phenomenology 
was carried out in Refs. [30, 34J, in which the implications of this Yukawa texture for the charged Higgs boson 
properties (masses and couplings), as well as the resulting pattern of charged Higgs boson decays and main production 
mechanisms at the LHC are discussed. Here, we will focus on Higgs boson pair production at 77 colliders through 
the 77 — >■ (f)i(t)j reaction,^ within the context of this more general version of the two-Higgs doublet model (2HDM-III). 
We will study the three distinct modes for Higgs boson pair production that are allowed: 77 — > AA, 77 — >■ A(j)a, and 
77 — >■ 0a0h, with (j)a = h,H. The impact of the Higgs potential on the production of pair s of light Higgs bosons, 
through the 77 hh reaction, has been studied recently [H, HI, H^l- The work of Ref. [3^ is focused to study the 
particular scenario where h couples to gauge bosons and fermions as the SM Higgs particle does, and finds that the 
cross section for the 77 — hh process can be much larger than the SM prediction. The same scenario was analyzed 
in Ref. [3^ . where it is taken into account the one-loop correction to the hhh vertex. It was found that this effect, 
together with the one-loop charged Higgs boson contribution, produces a considerable enhancement of the cross 
section. On the other hand, resonant effects due to the charged Higgs boson and the heavy neutral one H on the 
77 — 7> hh and 77 — ?> AA processes were studied in In this paper, we discuss the implications of assuming a 

four-zero Yukawa texture and general Higgs potential in the framework of the 2HDM-III on all the possible processes 
77 — > which include some scenarios of experimental interest did not considered in Refs. [H, HI, [131 ■ Apart 

from presenting analytical expressions for the three distinct reactions, namely 77 — > (j)a4'b,4'aA, AA^ we will analyze 
their unpolarized cross-sections in some scenarios that are still consistent with the most recent bounds on the model 
parameters as obtained from electroweak precision measurements [sl]. Higgs boson pair production at 77 colliders 
has already been studied in the SM by means of a nonlinear i?j-gauge j39|. and in the MSSM via a linear gauge 
[40| . Partial results were also obtained within the context of the 2HDM-II |4lj . 

On the other hand, the calculation of 77 — > scattering is far from trivial, mainly due to the plethora of Feynman 



1 The type-X (type-Y) 2HDM is referred as the type- IV (type-Ill) 2HDM in Ref. [H and as the type-P (type-IP) 2HDM in Ref. [13 . [l ^ . 
Sometimes the most general 2HDM, in which each fermion couples to both Higgs doublet fields, is called the type III 2HDM [17|| . All 
variants of the 2HDM were discussed recently in [Tsjl. The version that will be considered in this work will be called simply the 2HDM-III 
model. 

^ Throughout the paper, the symbol i}>i will stand for any of the three neutral Higgs bosons of the 2HDM, h, H, A, whereas flia will denote 
exclusively a CP-even Higgs boson. 
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diagrams arising in the gauge sector. Such a hard task can be significantly reheved if one uses an ap prop riate gauge- 
fixing procedure for the W boson. Rather than using the conventional linear gauge-fixing procedure [42] , we will use 
a nonlinear scheme that is covariant under the electromagnetic group f/e(l) [43]. We will use a general renormalizable 
nonlinear gauge-fixing procedure for the 2HDM fiS], which is intended to remove the most unphysical vertices from 
the interaction Lagrangian, thereby facilitating the calculation of radiative corrections considerably. It turns out that 
this class of gauges is tailored for the calculation of 77 — > <j>i(l)j scattering. We will show that such a gauge not only 
reduces dramatically the number of Feynman diagrams, but renders manifestly gauge-invariant and ultraviolet-finite 
amplitudes. The relevance of nonlinear gauges ^43)] to the calculation of radiative corrections has been emphasized 
by several authors not only within the context of the SM jl^, H^, |4^, but in some of its extensions such as the 
2HDM [13], the so-called 331 model [3, and also in the model independent effective Lagrangian approach 

The paper is organized as follows. In Sec. [TTjthe main features of the Higgs- Yukawa sector of the 2IIDM-III model 
are discussed. In Sec. IIIIl the main ingredients of the nonlinear i?j-gauge appropriate for this model are presented. 
Sec. IIVI is devoted to the cross-sections for 77 (j)i(f>j scattering, whereas numerical results are analyzed in Sec. [Vj 
Finally, the conclusions are presented in Sec. I VII 



II. THE HIGGS- YUKAWA SECTOR OF THE GENERAL TWO HIGSS DOUBLET MODEL 

In this section, we will discuss the main features of the general Higgs potential and the implementation of a specific 
four-zero texture in the Yukawa matrices within the 2IIDM-III. When a flavor symmetry in the Yukawa sector is 
implemented, discrete symmetries in the Higgs potential are not needed, so the most general Higgs potential must be 
introduced. 



A. The general Higgs potential in the 2HDM-III 

The 2HDM incorporates two scalar doublets of hypercharge -1-1: = and $2 — {(1)2 ,(f>2*)- The most 

general gauge invariant and CP-conserving potential can be written as [13] 

V{^i,^2) = M?!*!^'!) + M2(*2*2) - (/i?2(«'I*2) + H.c.) + ^X,{^\cS>,)^ + lAal^^^a)' + A3($t$i)($^ $2) 

+ A4($I$2)(*5$l) + Qa5($I$2)' + (A6($I$i) + A7($^$2)) (*I*2) + H.C. j , (1) 

where all parameters are assumed to be real^. In many discussions of the 2HDM, the terms proportional to Ag and A7 
are absent, as happens in the 2HDM type I and II where the discrete symmetry <I>i — ?> $1 and $2 — >■ —^2 is imposed in 
order to avoid dangerous flavor changing neutral current (FCNC) effects. However, in our model where mass matrices 
with a four-texture are considered, it is not necessary to implement the above discrete symmetry. Therefore, we must 
keep the terms proportional to Xq and A7. As we will show below, these parameters play an important role in the 
77 — 4>i4>j reactions. It is worth commenting that the parameters Ag and A7 are essential to obtain the decoupling 
limit of the model in which only one CP-even scalar is light. As long as these terms exist, there are two independent 
energy scales, v and Athdm, and the spectrum of Higgs boson masses is such that rrih is of the order of u, whereas 
niH, iTiA and m//± are all of the order of Athdm [13] • In this case, all of the heavy Higgs bosons decouple in the limit 
of Athdm ^ v, according to the decoupling theorem. On the other hand, when the scalar potential does respect the 
discrete symmetry, it is impossible to have two independent energy scales [13] • As a consequence, all of the physical 
scalar masses lie on the Fermi scale v. Since v is already fixed by the experiment, a very heavy Higgs boson can only 
arise through a large dimensionless coupling constant A^. In this scenario the decoupling theorem is no longer valid, 
thereby opening the possibility for the appearance of nondecoupling effects. In addition, since the scalar potential 
contains some terms that violate the SU (2) custodial symmetry, nondecoupling effects can arise in one-loop induced 
Higgs boson couplings 'si']. 

The scalar potential (|T]) has been diagonalized to yield the mass-eigenstates fields. The charged components of the 
doublets lead to a physical charged Higgs boson and the pseudo-Goldstone boson associated with the W gauge field: 

Gi, = ^tcp+4sp, (2) 
= -cbfsp+cbfcp, (3) 



The X(i and A7 parameters are complex in general, but we will assume that they are real by simplicity. 
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with tan/3 = V2/V1 = tp, being vi/\/2 (W2/V2) the vacuum expectation vakie (VEV) associated with $1 ($2), and 

„2 2 

^ - -v^{Xi + X5+t/Xe + tpXj), (4) 

S/3C/3 2 

where we have introduced the shorthand notation = sin/3 and = cos/3. On the other hand, the imaginary part 
of the neutral components (p'^j defines the neutral CP-odd scalar and the pseudo-Goldstone boson associated with the 
Z gauge boson. The corresponding rotation is given by 



Gz = 0?/C^ + 0§/S/3, (5) 
A^-(PljSp + 4'2iC^, (6) 



where 



m\ = mjj± +^v^{X4- X5). (7) 

Finally, the real part of the neutral components of the (f)^j^ doublets defines the CP-even Higgs bosons h and H. The 
mass matrix has the form: 

Mn. = ( "^'^ ] , (8) 

where 

mil = rn^sp + w^(Aic^ + s^X^ + 2spcpXG), (9) 

^22 = rn\cl + w^(A2s| -|- cjX^ + 2spcpX7), (10) 

mi2 = -rri^spcf) + {{Xj, + A4).s^c^ + Aec^ + X-jsj^ . (11) 
The physical CP-even states, h and H , are written as 

^ = Ac«+<^2i?Sa, (12) 

= +02flCa, (13) 

where 



tan 2a 



2mi2 

TOii - m22 ' 



(14) 



and 



TO22 ± \ [mil - "^22) 



4m22 



(15) 



B. The Yukawa sector in the 2HDM-III with a four-zero texture 

We shall follow Refs. [1^ H^, where a specific four-zero texture has been implemented for the Yukawa matrices 
within the 2HDM-III. This allows one to express the couplings of the neutral and charged Higgs bosons in terms of 
the fermion masses, Cabibbo-Kobayashi-Maskawa (CKM) mixing angles and certain dimensionless parameters, which 
are to be bounded by current experimental constraints. The Yukawa Lagrangian is written as follows: 

Cy = Y^Ql^iUr + Y^Ql^2Ur + YfQL^idR + Y2^QL^2dR + H.c, (16) 

where $1,2 = i4'i2^ ^'1,2)'^ refer to the two Higgs doublets, $1,2 = ^(^2^1.2^ Ql denotes the left-handed quark doublet, 

ur and dR are the right-handed quarks singlets, and Y^^ denotes the (3 x 3) Yukawa matrices. The Yukawa lepton 
sector is given by a similar expression. 

Since the fermionic contribution to the one-loop reactions 77 — > 4'i4'j is given by vertices ff4>i involving 
only neutral Higgs bosons, we will concentrate only in this parte of the Yukawa sector. After implementing the 
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diagonalizations carried out in the Higgs potential and in the Yukawa sector^, the interactions of the neutral Higgs 
bosons {h'^ , , A^) with quark pairs acquire the following form; 



'-Y 



9 f ITT'd, 



2 \mw 
, 9 f 



cos a 



\/2 sin(c 



( ^ 1 {yA 



2 \mw 

W ( rrid, 



2 \mw 

9 f ruui 
2 \mw 

9 ( 



sma 

a°i3 

COS p 



tan/3(5i 



g COS P \ rud. 

\/2 cos(a ~ (i) f rriw 
g cos f3 \ rudi 

\/2 / mw 



2 J^J 



g cos 13 \ iTLdi 



2 JiJ 



sina V2 sinja - f3) ( mw , 



2 \mw 
^9 ( 



cos a 
sm p 



cot/3(5i. 



g sin /3 \ r7i„ 

-\/2 cos(q! — (3) ( niw 
g sin /3 



5sin/3U«J^ ' 
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(17) 



where i = 1, 2, 3, with di ^ d,d2 = s, c?3 = b,ui = u, U2 — c,U3 = t. These couplings depend on the rotated matrices 
— Vq [n — 1 when q = u, and n ~ 2 when q — d ). Here Vq is the diagonalizing mass matrix. In order to 

evaluate Y^ we need to focus in the quark mass matrix, which is given by, 



Ml 



V2 



{v^Y,' +V2Y^^), {q = u,d). 



(18) 



We shall consider that all Yukawa matrices have the Hermitian four-zero texture form [si'l, and the quark masses 
have the same form, which are given by: 



Cq 
M" = I C* Bq Bq 

B; Aq 



{q = u,d). 



(19) 



This is called a four-zero texture because one assumes that the Yukawa matrices are Hermitian, therefore each u and 
d type Yukawa matrix contains two independent zeros. According to current analysis this type of texture satisfies 
the experimental constraints (i.e. the Flavor Violating Higgs interaction) and at the same time it permits to derive 
analytical expressions for the Higgs boson fermion couplings [2^, [H, Is^, [sl, [s^ . 



To diagonalize these mass matrices, we use the matrix Vq in the following way: 

M9 = VqWV^. 



(20) 



Following the analysis in one can derive a better approximation for the product Vq Y^ V^ , expressing the rotated 



matrix Y^, in the form 



mfm'j 



[X 



01] 



mfrrij 



(21) 



where x's are unknown dimensionless parameters of the model, they come from the election of a specific texture of the 
Yukawa matrices. In order to perform our phenomenological study, we find it convenient to rewrite the Lagrangian 



* The details of both diagonalizations are presented in Ref. [Ssll . 

^ Vq is built as a product of two matrices Oq and Pq in the form Vq = Pq . These matrices are given in [3l| . 
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given in Eq. p7|) in terms of the coefficients [Xnli,-; as follows: 



md. 



cos a sm{a ~ f3) / ^JmJrrQ,. > 
TOvv ) cosp 



Tow 



*5 J 
— d, 

2 



: Ml 



5 - 
2 

«5 - 

2 



\ TOvK / cos /3 



cos/3 

sin a cos(a - /3) / v^WdlTO^X ^ 

^ 72 cos /3 I mvi. j 
1 / y/frTdjrQ 



m„. \ sin a 



tan /3 (5i 



TOvy / sin /3 



cos V '^VF 
sin(a - 13) f ./n^^m, 



Xij 



djh° 



V2 sin/3 



mw 



TO,,, \ cos a cos(a ~ P) f -/^^W^ti^ 



TOvK / sm 



cot /3 (5y 



\/2 sin /3 V "^w 

1 / ^TO„,TO„^ 



V2 sin/? 



TOVF 



Xi 



(22) 



where we have redefined [xi]^ = x^j 



and [x2Jij = Xij- it was discussed in Ref. [28], most low-energy processes 
imply weak bounds on the coefficients xfj, which turn out to be of 0(1). Based on the analysis oi B ^ Xs^ [52|.|53|. 

we find the bounds: Ixsa'^l ^ 1 for 0.1 < tan/3 < 70 [s^. Other constraints on the charged Higgs mass and tan/3, can 
be obtained from anomalous magnetic moment of the muon Aa^, the p parameter, as well as B-decays into the tau 
lepton [sj,!!^. For instance, as it can be read from Ref. [s^, one has that the decay B tv, implies a constraint 
such that for m/^+ ~ 200 (300) GeV, values of tan/3 less than about 30 (50) are still allowed, within MSSM or 
THDM-II. However, these constraints can only be taken as estimates, as it is likely that they would be modified for 
THDM-III. A more detailed analysis that includes the most recent data is underway [13, [11]. On the other hand. 



the condition 



- H + 



< 



in the frame of the 2HDM-II implies 



However, we found that in 2HDM-HI 



- ff + 



'-H + 
™ff+ 

1 



4V2ir tan /32 



which leads to 0.3 < tan/3 < 130. 



[30|, we have checked numerically that this 



leads to 0.08 < tan/3 < 200 when jxasJ ~ 1 and 0.3 < tan/3 < 130 as long as jxasJ recovering the result for the 
case of the 2HDM-n [H [Hi . 

Other important bounds on Jxsal and tan/? come from radiative corrections to the process T(Z —> 66), specially 
the hadronic branching fraction of Z bosons to 66 (Rb) and the 6 quark asymmetry (At) impossed a high restriction 
[13, [iO]- We can get bounds for tan/3: in the case X33 = 1 and mH+ ~ 200(300) GeV, the range tan/3 > 0.3(0.2) is 
allowed, while in the scenario X33'' = — 1 and to//+ ~ 200(300) GeV, tan^ > 5(3) is permitted. 

On the other hand, the leading contribution to Bq — Bq mixing in the regime small tan /3 is given by the charged 
Higgs sector. Following the Ref. ^[ll, we get for the case X33'' = 1 and mH+ 200(300) GeV, tan/3 > 0.2(0.25) 
is allowed. Combining the criteria of the analysis for the radiative corrections of Zbb vertex and Bq — Bq mixing, 
tan/? > 0.3 is allowed for > 170 GcV and Xss^ — 1- However, when = ^1 and 111^+ < 600 GeV, tan/? < 2 

is disfavored. 

Besides, following the analysis of the Ref. [H, [H, [llj, one can get the deviation Apo of the parameter po — 
/ pM'^C'^ of our version 2HDM-HI, where the p in the denominator absorbs all the SM corrections, and the most 
important SM correction at 1-loop level comes from the heavy top-quark. According the reported value of po is [6]| 



1 nnA +0-0027 ,„ , 

Pa = 1.004 -0.0007 (2cr) . 
In terms of new physics (2HDM-HI) the constraint becomes: 

- 0.0007 < Ap2HDM-iii < 0.0027 . 
In 2HDM Po receives contribution from the Higgs bosons given by, in the context of model HI [l3, [Hi 

Gp 



(23) 



(24) 



Ap 



2HDM-III — 



sin'ia- P)F{Mh±, Ma, Mho) 



+ cos^ia- P)F{Mh±, Ma, Mho) 



(25) 
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where 



F {nil, 7112, ms) = 




Since pg is constrained to be around 1 we have to minimize the contributions of Ap2HDM-iii- This is obtained for 
the case a — 0,7r/2, and the parameter space of the scalar sector is strongly reduced when decoupling between Higgs 
bosons, i.e. Am^j — mi — rrij > 100 GeV (m^ = ni/jo, m^/o, to^o, m^±). However, is possible to avoid the constraint 
for Ap2HDM-iii if the decouphng source Aniij ^ 20 GeV or Am^ ~ 100 GeV and one Higgs very heavy (e.g. to/^o > 1 
TeV). When a = /3 ± 7r/2 the allowed parameter region is larger and one can avoid the constraints of the p parameter 
with or without decoupling. Another interesting possibility is when we have the case quasi-degenerate between the 
masses of CP-even Higgs boson (H) and the charged Higgs boson {H^). The reason is that, in a 2HDM, the custodial 
symmetry may be implemented either with rriA = mff± or with mn — rnH± [13] ■ The full study of the p parameter 
in our version 2HDM-III will be presented elsewhere [57|]. 

Hereafter, we shall refer to three benchmark scenarios, namely: 

• Scenario I (the decoupling limit). In this scenario, h assumes the role of the SM Higgs boson hsM and is 
essentially independent from the diverse versions of the model. We have chosen to discuss the 77 hh process 
within this context to illustrate the decoupling nature of the heavy Higgs effects [s^. We will take nih = 120 
GeV and toa ^ "i^i ^ uih >> v. Two cases will be considered, one when the parameter of the Higgs potential 
/X12 is of the order of the Fermi scale, /ii2 ~ v, and other when this parameter is much larger that such scale, 
H12 » V. 

• Scenario II ( SM-like) "sO*]. This is a scenario of the 2HDM-HI in which the couplings hVV {V = W,Z), 
hhh, hhhh are nearly indistinguishable from the corresponding hsM^ whereas the hf f couplings can deviate 
significantly from the corresponding hsAiff ones. We will take the following values for the parameters of the 
model mil = 120 GeV, = 110 GeV, fiu — 130 GeV, m/^± ^ mn ^ mA+mt GeV, and a — /3±7r/2. Within 
the two subscenarios: (A7 = — Ag = —0.1) and (A7 = — Ag = —1). 

• Scenario III (a more general case of 2HDM-III). In this scenario, more general couplings of neutral 
Higgs bosons to SM particles are assumed. We will include the contributions of the parameters of the Higgs 
potential Ag and A7, as well as the contributions of the Yukawa texture in the couplings </>//. Within this 
scenario, a degenerate case, a nondegenerate case and the case with a light CP-odd scalar will be considered. 
In the nondegenerate case, we will choose mu± — 400 GeV, mA = 350 GeV, mn — 520 GeV, mh — 120 
GeV, and p,i2 = 120 GeV. On the other side, in the degenerate case we will choose mH± = itih — ttla — 300 
GeV, with /ii2 = 60 GeV and m^ — 120 GeV. For the case with light CP-odd scalar, we will choose = 50 
GeV, mh = 120 GeV, m^i = 350 GeV, mj^ = 400 GeV, and H12 = 70 GeV. In ah cases, the a = /3 and 
a — (3 ± tt/2 possibilities will be considered. In the nondegenerate case, the set of values Ay = — Ag = — 1 and 
A7 — — Ag = —0.1 will be considered. As we will see below, only the values Ay = — Ag = —1 are relevant for the 
degenerate case and when we study the case with a light CP-odd scalar. 

In all the above scenarios, we consider the constraints imposed by perturbativity, Z — >■ bb, po parameter, and 
— mixing. Our predictions will be consistent with current bounds on the charged Higgs mass obtained at 
Tevatron [63] and LEP2 [6l|, 113] , as well as with those derived theoretically [65j . 

III. THE GAUGE-FIXING PROCEDURE 

As already mentioned, in calculating the 77 — >■ (l)i(f>j reaction we will define the W gauge boson propagator using 
a nonlinear gauge-fixing procedure that is covariant under the electromagnetic gauge group and consistent with 
renormalization theory. The details of this gauge-fixing procedure for the 2HDM has been reported recently in 0| . 
Here, we present the gauge-fixing functions, including some results and comments that are needed in calculating the 
amplitude for the 77 — process. 

To begin with, we discuss the most general structure of gauge-fixing functions /" and / for the SUl{2) and f/y(l) 
gauge groups that are allowed by renormalization theory and covariance under the electromagnetic gauge group in 
the context of the 2HDM. As it is stressed in Ref.fiS], our main aim is to remove the most nonphysical vertices that 
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are generated by the Higgs kinetic-energy term. This is achieved by introducing the following nonlinear gauge-fixing 
functions 

r = fv + fs, (26) 
/ - /y + fs, (27) 



where 



and 



= (<5'^'a^-.gV'^'i?M)w^'^ (28) 

= E ['^I(^" - ^^"''<'>o^ - + le'^^'a')^.- 

1=1 

+ie3"^(c/3$l + sp<i>l)a\cp<S>i + 5^$2)}, (29) 



fv = (30) 

i=l 

In the above expressions, = (0,Ui/A/2), cr° are the Pauli matrices, and and i?^ are the gauge fields associated 
with the electroweak group. Our gauge-fixing functions contain the conventional linear functions as a particular case, 
which are obtained when e^"^ is set to zero. Also, it is worth mentioning that this gauge-fixing procedure contains as 
a particular case an analogous gauge scheme for the minimal SM jiB], which becomes evident when the $i doublet is 
associated with the SM one and (3 is set to zero. 

To fully appreciate the structure of the gauge-fixing functions, it is convenient to express them in terms of mass 
eigenstates fields. After this, one obtains for the vector sector 

(32) 
(33) 
(34) 

and for the scalar sector 

/5 =-^(^"-*GzjG+, (35) 

/I = -CmzGz, (36) 
/# = 0, (37) 

where ip^ = v + cp^aH + sp^^h and — d^ — ig'B^, being g' the coupling constant associated with the ?7y(l) group. 
We can see that both fy and fg are nonlinear and transform covariantly under the ?7e(l) group, as contains the 
covariant derivative associated with this group. 

Following the study of Ref. [Si, the gauge- fixing Lagrangian, £_b, can then be written as 



n = 




= 








2 





where 



r 

^-BSV — 



£b — ^BV + ^BS + ^BSV, (38) 

cbv = ~^{D,w+nHD.w+n - Y^id.zn' ^{d,An\ m 

Cbs - + ^^^"^ - ^Cm|G|, (40) 

{D^W+^)H^^ - iGz)G+, - {D^W+>'){v + iGz)G^] + mzGzd^Z^'. (41) 
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We restrict our discussion to present some comments concerning the impact of the Cbv ^ ^bS: a-nd Cbsv Lagrangians 
on the Yang-Mills, the Higgs kinetic-energy, and the Higgs potential sectors, respectively. First of all, the term Cbv 
defines the propagators of the gauge fields and also modifies nontrivially the Lorentz structure of the trilinear and 
quartic vertices arising from the Yang-Mills sector. Indeed, with the exception of the WWWW vertex, all trilinear 
and quartic vertices are modified by the gauge fixing procedure. Since the term that introduces the modifications 
on these vertices is invariant under the Ue{l) group [4J], the trilinear electromagnetic vertices satisfy QED-like Ward 
identities. This fact is relevant for radiative corrections as such a symmetry greatly simplifies this class of calculations. 

As for the Cbs term, it defines the masses of the Gw and Gz fields and modifies some nonphysical couplings 
arising from the Higgs potential. In this gauge, the couplings between scalar fields arise solely from the sum of the 
following terms —V{^i, ^2) + ^bs- The last term in this sum leads to modifications in the strength of the nonphysical 
couplings iJG^G^, hG^G^, H^G^G^, K^G^G^, HhG^G^, and G%G^G^. The physical couplings remain 
unchanged, as required. 

On the other hand, the term Cbsv considerably affects the Higgs kinetic-energy sector of the theory since it removes 
several nonphysical vertices. When these two terms are combined, one finds |44l | that not only the the mixing terms 
W — Gw and Z — Gz are removed from the theory, as it occurs in conventional linear gauges, but also the nonphysical 
vertices WGw7, WGwZ, HWGwl, hWGwl, GzWGwl, HWGwZ, hWGwZ, and GzWGwZ are removed. In 
addition, the unphysical vertices HWGw, hWGw, and GzWGw are modified. Once again, it should be emphasized 
that the couplings involving only physical scalars are not modified by the gauge-fixing procedure. 

As far as the ghost sector is concerned, it is shown in Ref. [44] that it shows new interesting aspects that are not 
present in conventional linear gauges, such as manifest electromagnetic gauge invariance and the presence of quartic 
ghost interactions. As a consequence of [/e(l)-gauge invariance, the corresponding electromagnetic couplings satisfy 
QED-like Ward identities, which considerably simplifies the loop calculations associated with the 77 — 0^0^ process. 

The gauge-dependent Feynman rules necessary for the calculation of the 77 — > (piipj processes are not presented 
here, as they are given in Ref. [44|]. The rest of Feynman rules, which do not depend of the gauge-fixing procedure, 
are given in an Appendix |^ 



IV. THE PROCESSES 77 0,0^ 

We now will exploit the nonlinear i?^-gauge already introduced to calculate the 77 — (j}i4ij process. To begin 
with, we would like to discuss the basics of this process, such as its kinematics and the gauge structure dictated by 
electromagnetic gauge invariance. To this end, we use the following notation: 

A^(fci) + A,{k2) Mh) + (t>j{k4), (42) 



where the particle momenta satisfy the kinematic relation ki + k2 — + k^. The Mandelstam variables associated 

'4+ki 



with this process are s = {ki + ^2)^, t = (^2 — ^3)^, and u = (fci — k^)"^, which fulfill the relationship s + t + u = k'^ + k'l 



One useful quantity is the transversal momentum, given by 



k^^-{tu-klkl). (43) 



s 



As far as the gauge structures are concerned, two possibilities arise depending on the CP properties of the final 
particles. One possibility corresponds to final states with two Higgs particles both CP-even, i.e. (l>a<Pa or both CP- 
odd, AA. Although it is possible to construct at least three gauge electromagnetic structures, only two of them are 
independent. We find it convenient to use the following basis: 



Pl^Ll' = — (k2f^ki„ ~ ki ■ k2gi_,v), (44) 
^/pS - 2 



-P2^iy = T2— -;^krpsgi_,„ + k^k2f,ki^ - 2^2 • fcafca^fci,, - 2fci • fc3fc2^fc3^ -I- 2fci • k2k3f^k3^ , (45) 



which is orthonormal in the sense that P^'^ Pji^^^ — 5ij. Another possibility corresponds to final particles with distinct 
CP properties, i.e., Acjja- The corresponding gauge structures can be assembled by combining the Levi-Civita tensor 
and the 4- vectors fci^, k2^, and /ca^. However, not all the combinations are independent as some of them can be 
eliminated with the help of Schouten's identity. Once those redundant structures are removed, we are left with two 
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independent gauge structures. We choose the foUowing basis 



Plfiu = —<^ii,^al3kikt^, (46) 



2V2 



2^1^ — 1.2 



krp S 



which also are orthonormal. The invariant amplitude can be written as follows: 

M=M^,e''iki,Xiy{k2,X2), (48) 
where e'^(A:i, Ai) and e'^(A:2, A2) are the polarization vectors of the photons. The tensor amplitude reads 



■^5 V^a Y-'a ; 



AiPif,^ + A2P2f,^ for AA, 

AlPi^^ + ^2^2^-^ for A(l>a 



V2s 



w 



and the unpolarized cross-section for the process 77 — ^ (t>i4>i is given by 

rto 



spin 

i {\A,\' + \A2n 



USns^s^el ~ ~ ^^^^ 



where e = 2 if the final particles are identical and 1 otherwise. The integration limits are 



(51) 



We will present below the amplitudes for the three available processes: 77 — >■ AA, 77 — > A(l)a, and 77 — >■ (j)a(t>b- The 
absence of the unphysical jWGw and ipajWGw vertices ^ introduces considerable simplifications in the calculations. 
In particular, there are a significant reduction in the number of diagrams with respect to those appearing in a linear 
gauge. Also, the contributions can be grouped into distinct sets of diagrams which lead to finite and gauge invariant 
results by their own (see Figs. [T][5]). 

We find it convenient to express our results in terms of Passarino-Veltman scalar functions. For those scalar 
functions arising from the sets of diagrams of Figs. [1] through |4l an unambiguous shorthand notation can be used: 

Co{a,b) = Coikl,klika + kbf,m^,m^,m^), (52) 
Do(a, 6, c) = DQ{kl,kl,kl,{ka + kb + kcf,{ka + khf,{kb + kcY,m^,m'^,m^,m^), (53) 

where the right-hand side has been expressed in the notation of Ref. [g^. In addition, a, b, and c run over 1 to 3, 
and m is the mass of the particle circulating in the loop. Below it will be evident why it is not necessary to specify 
m as an argument of the scalar functions. Unfortunately this scheme cannot be used for the scalar functions arising 
from Fig. [5] since two different particles circulate in the loop, so an adequate notation for these scalar functions will 
be given below. 

It is also convenient to introduce the following dimensionless variables Xf = m'j/s, xa = rr^^js, Xa = Tn^js {iria 
stands for the mass of (pa), xh = fn^i/s, xw — m'w/s, xt = k'^/s, Xt — t/s, Xu — u/s, and 7^ — Ti/^/s, being 
the total width of (pi . 

A. The process 77 — > AA 

This process receives contributions from charged fermions, the W gauge boson and the charged Higgs boson H^. 
These contributions arise through the diagrams shown in Figs. [U [21 [31 and [H There are not contributions from 



The A-yWGw vertex is not generated by the theory. 
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FIG. 1: Feynman diagrams contributing to the 77 — >■ 4'i4'j processes in the nonlinear gauge, in which this contribution is 
ultraviolet finite and gauge invariant by itself. In the case of the W contribution, similar diagrams are to be considered for the 
pseudo-Goldstone and ghost fields. Crossed diagrams are not shown. 




(i) (ii) 
FIG. 2: The same as in Fig. [T] 



diagrams of Fig. |4] since the AWGw vertex is not generated by the theory. As pointed out before, in the nonlinear 
gauge, each set of diagrams leads to an ultraviolet finite and gauge invariant amplitude. All of the diagrams of Figs. 
[U [21 [3l and [5] contribute to the Ai amplitude, but only the diagrams of Figs. |3] and [5] contribute to A2. These 
contributions can conveniently written as: 



Ai = All 
Ai = A23 



A12 

-425, 



^13+^15, 



(54) 
(55) 



where the subscript j in Aij denotes the contributions coming from each set of diagrams, which can be expressed in 
terms of the contributions of each kind of particles circulating in the loops as follows: 



All = Fii + Sii + Gil, 

A12 = S12 + Gi2, 
^13 — Fi'i, 

A23 ~ F23, 



(56) 
(57) 
(58) 
(59) 



where we will use the letters Fij, Sij, and Gij to denote the contributions coming from fermion, scalar, and gauge 
particles, respectively. The set of diagrams in Fig. [T] yield 



Fii 
Sii ■ 
Gil 



f 4>a 



0a 



1- Xa 
Xw 



-[2 + (4a;/-l),sCo(l,2)], 



GtpaWwGci)^ 



AA 



l-Xa+ iyfxl'^a 

XW_ 

1 - + iJx^^a 



[l+2XffSCo(l,2)], 



12 + 



xw 



(l22;vy -4 + a;a)2sCo(l,2) 



(60) 
(61) 
(62) 



^AAAAA^ 



^AAAAA> 



f, W^H^ 




(11) 



FIG. 3: The same as in Fig. [T] 
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FIG. 4: Feynman diagrams contributing to the 77 — >■ (j>i4'j processes in the nonhnear gauge. The contribution of these diagrams 
is finite and gauge invariant by itself. Crossed diagrams are not shown. 



where Nf is the color index and Q/ is the electric charge in units of the positron charge, Q,p^(f f.vv,<p<j>) ^-re functions of 
the couplings g<t>a(f f,vv,4>diXi,^^^'^^ they are given in Appendix. Notice that there is no contribution from the one-loop 
77Z* off-shell coupling [67| as the Z boson only couples to pairs of the form Acjja ■ In the set of diagrams of Fig. [2] 
there are only contributions from the charged scalar boson and the W gauge boson. The corresponding amplitudes 



are 



5^, = _(^^™A)[i + 2xhsCo(1,2)], 



G 



12 



1 - (1 - 2xw)sCo{l, 2) + ( ^"'g^'" ) [I + 2xwsCo{l, 2)] 



(63) 
(64) 



There are contributions to G12 from the W boson and its associated pseudo-Goldstone boson, but not from the ghost 
field because of the absence of the C^C^AA vertex. As far as the diagrams of Fig. |3] is concerned, there are no 
contributions from scalar or gauge particles due to the absence of the H^H^A and WW A vertices. The fermion 
contributions are given by 



i^i3 = 2 ^ NfQ}glj^ (^) \ 2 + 4a;/sCo(l, 2) + 2xas \{xt - xa)Co{2, 3) + (x„ - xa)Co{1, 3) 

J \ X]^ / y L 



xaxtDo{1, 3, 2) + xf{l - 2xa) -Do(l, 2, 3) + D„{1, 3, 2) + Do{2, 1, 3) 



+ 



(65) 



23 



1 - 2XA 



2XT 



s 

2^ 



Coil, 2) + i2xA - l)Co(3, 4) + (xt - XA)Coi2, 3) + (x„ - XA)Coil, 3) 
2xt{xt - xa)Co{2,3) + 2xu{xu - a;^)Co(l,3) - (1 - 2x^)^0(1, 2) - {xt - x„)2Co(3,4) 
2(1 - 2xA)xfXT \Do{1, 2, 3) + Do{l, 3, 2) + Do{2, 1, 3) 



xtix^t + x\)Do{l, 2, 3) + xjxl + x\)Do{2, 1, 3) 



(66) 
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vwwv 

VW\A/\> 



WWW 

'ww\A^ 



www- 
^A/wv\^- 



' G,, 




vwwv 

H 



-\A/WW 





FIG. 5: The same as in Fig. H 



Finally, the contributions from Fig. [5] read 

12 6 



a— 3 



A 



25 



12 



(67) 
(68) 
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where 



fl 



9\ 



9\ 



95 



-8(1 - xw), 



9l 
9l 
9l 



fll = fl2 = ^o(a;n -Xa), 



fl = /s = /g = /lo = M^t - xa), 

Xw 



{2xw + Ao), 



\xw ' 

2{2xH+Ao), 

{xT + XH + Xw)Ao 



-Ax 



H 



Xw 



(69) 
(70) 

(71) 

(72) 

(73) 
(74) 
(75) 



fl = 
fl - 



fl 
fl 
fl 



/I = {2xw + ^o) 2xT + (a;* - Xuf , 

{2xw + Ao){2xw - 2xh - Xt - Xu) 

{2xw + Ao){2xh — 2xw — Xt — Xy) — SxtXw] 



fl 



fll = fi2 = (xA - Xu) XuAq- 2xw{xa ~ XuY 



— fs — fg ^ fio ^ fb ^ ^t)' 



gl = (Ao + 2xw)[2xtXh + {xh +Xu- xw) 



(76) 

(77) 
(78) 

(79) 

(80) 

(81) 



gl = Axw{xa - Xu)^ - Axw{xa - Xuf + {xa - XuT 2(2a:^ + 2xh - ^xw + '^)xw + Aq{1 - 2xw) 
+2{xA - Xu) 2xw{xw - xa) + Aq{xh - xa) + + "^xwixA + 1)^0 + '^xaxhxw, 



(82) 



with 



53 = 9i{xu Xt), 

gl = 52(2;* ^ Xu), 



2 2 

95 = 9e = 



i|(Ao + 2xw){2xh ~ 2xw + xt+ XuT 

+2[Ao(a;T - xt- Xu) - 2xw{xt + Xu)] {2xh - 2xw + xt + Xu) 
+8xtxw{Aq + 2xw) + Sx^xw — 2Ao{xt + Xu)xt 
+ {Ao + 2xw){xt+Xu)^y 



Aq ^ x\- 2{XH + xw)xA + {xH - xw)'^, 
The arguments for the Co (a) and D(){a) scalar functfons are presented m Appendix B. 



(83) 
(84) 



(85) 



(86) 



B. The process 77 — > 0a 06 

We turn now to the amphtudes associated with the final states hh, hH, and HH. There are contributions from all 
the sets of diagrams shown in Figs. [T]to[5] According to our notation, the contributions can be organized as follows: 



Ai = An + A12 + Ai3 + An + Ai5, 

^2 =^23+^24+^25, 



(87) 
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where 



All = Fii + Sii + Gil, 

A-12 ~ Si2 + Gi2, 

Ai3 — Fi3 + 5i3 + Gi3, 

A-23 — F23 + S23 + G23. 



(89) 
(90) 
(91) 

(92) 



The Ai4, A15 A24, and A25 partial amphtudes only receive contributions from the pairs {W,Gw) and {W,H^). The 
contributions coming from the set of diagrams in Fig. [T]are given by 



Fii = Y.NfQ}Y,{g,j^g^^,,^j- 

f 4>c 

^"--z^i i 

G"--Z.i 2 



Xw 



[2 + (4x/-l)sCo(l,2)], 



1- Xc + i^cli 

Xw 



1 - + iy/x^lc 



[l + 2a;ifsCo(l,2)], 



12 + — + 2(-4 + 12a;w + Xc)sGo{l, 2) 

Xw 



(93) 
(94) 
(95) 



^12 = -^Il±Ipt^[l + 2xHsC,{l, 2)], 



G12 = 8sCo(l, 2)5'"' + (2Co(l, 2)sxw + 1)(2^0 



tWW - Tk'G* a 



sab\ 



(96) 
(97) 



^13 



E^/^?/(^0./7^0./7)(^)|4[l + 2^/-^^o(l,2) 



(1 + + x„ - 8a;/)s ^(a;„ - Xa)Co{l, 3) + (a;„ - Xb)Co{2, 4) + (.t* - Xa)Go{2, 3) + (a;* - Xb)Go{l, 4) 
2x/(l + Xa + a;b - 8xf)s^ [i^o(l, 2, 3) + Do{2, 1, 3) 



16xj + 2a;/(4a;T 



Xb-l)-XT{Xa+Xb) s^Do(l,3,2) ^, 



(98) 



'13 



^xwsl [xa - x„)Co(l,3) + {xb - a::„)Co(2,4) + {xa - a;t)Co(2,3) + 



[Xb - Xt)Co(l, 4) + 2xHS i5o(l, 2, 3) + Do(2, 1,3) + Do(l, 3, 2) + xts^?o(1, 3, 2) 



(99) 
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bWW 



Xw 



siAc 



)C^o(2,4) + ( Xa Xf )Co(2,3) 



+ {xb - xt)Go{l, 4) + 2x1^5(^0 - ^xw) ^0(1, 2, 3) + -Do(2, 1, 3) + 7^o(l, 3, 2) 



sxTAGDo{l,i,2) 



+ 



(100) 



with 



Ag = 16a;^ + XbXa + 2a;vF(a;a + Xf,). 



(101) 



(40(^„ivw^0,^(,w^w')s| - 2Co(l,2) + 

{xa + Xb) [{xa - Xu)Cq{1, 3) + (xfe - a;„)Co(2, 4) + (a;„ - Xt)Co{2, 3) + {xb - Xt)Co{l, 4) 

2xw{xt + Xu)s \Do{1, 2, 3) + Do{2, 1, 3) + Do{l, 3, 2)1 + xrixa + Xb)sDo{l, 3, 2) L 



IZ D 

^15 = (Gw±ht4>^ 



a=3 



a=l 



fl = -AXH, 

fl = %{xw - 1), 

fr = fl = -{xa -xt)f, 
/g = /lo = (xb -xt)f, 

/ll = /12 = {Xb — Xu)f, 



9I =9l = 

92=94 = 



95 =96 



-2xH{f + 2), 
-2xw(f + 2^), 

\ Xw ' 

- (xh +Xw+XT^f - AxhXw , 



f = 



{xh - XwY - (xh + Xw)iXa + Xb) + XaXb 



Xw 



F23 = Y.^fQ'ASi>Jf)^ 
f 



XT 



Co (1,2) 



( 



x',+xl 2x,Xb \ ^ _ 8xf)Co{3, 4) + 

XT ' 
Xtixt - 8xf ) + XaXb' 



XT 



Xu{xu - 8a;/) + 



XT 



{xa - xt)Co{2, 3) + {xb - xt)Co{l, 4) + 



(xa - a;«)Co(l, 3) + (xb ~ x„)Co(2, 4) 



2xf{xt+Xu - 8x/)sDo(l,3,2) + 



16x} + 2xf{xt +Xu- 4-^) + {xl + XaXb)— Do{l, 2, 3) + 

■' XT XT-' 

16x} + 2xf{xt +Xu- 4^) + {xl + XaXb)—\ Dq{2, 1, 3) I, 

Xt Xt^ 



G23 = {G4,aWwG4>bWw)( ^Acs(a + xwB 

\XtX-w ^ 
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with A and B given by 

A = -{xt+ Xu)Co{l, 2) + (a;? +xl- 2xaXb)Co{3, 4) + xt{xa - Xt)Co{2, 3) + Xt{xb - Xt)Co{l, 4) + 
Xuixa - Xu)Cq{1, 3) + Xuixb - Xu)Co{2, 4) + s \x'^Dq{1, 2, 3) + a;2£)o(2, 1, 3) 



B = 2xts 



£>o(l,2,3) + £)o(2,l,3) + L»o(l,3,2) 



(117) 
(118) 



Finally, the A24 and ^25 amplitudes read 

+ 4 + 2a;aa;fe)Co(l, 2) - (x^ + a;^ _ 2xaa;6)(xt + a;„)Co(3, 4) 



(a;? + a;aa;6) (a;,, - a;t)Co(2, 3) + (a;^ - a;t)Co(l, 4) - {xl + XaXb) (a;„ - a;„)Co(l, 3) + {xb - a;„)Co(2, 4) 



S\ XT 



2xw{xt + Xu) + XT -Do(l, 3, 2) + - 2xwxt{xu + Xt) + a;| - Xt{xaXb + x^) -Do(l, 2, 3) 



2xwXt{xu + Xt) + x'^- Xu{xaXb + xl) Do{2, 1, 3) 



+ 



(119) 



where 



^25 = (pW±HT,j,aQw±HT^^^ (^) 
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X:/„'C/o(a) + s^5^i5o(a) 



a=l 



a=l 



/i = /I = (/ + 2) [(a;* - a;„)2 + 2a;T] 

fi = -U ^2){2xH-2xw+Xu^xt), 

!l = (/ + 2)(2.T// - 2x14^ - a;„ - a;*) - 8xT 

/b = /e ^ (-^'a ^ 2;„) [/a;^ 2(a:£i a;„)j , 

/? = /I = (2^1 ~ Xt) \!xt - 2{xa - xt){xb - Xt)] , 

/g = /lo = (a^b - Xt) [fxt - 2{xa - Xt){xb - Xt)] , 

/ll = /12 ~ {Xb ~ Xu) [fXu — 2[Xa — Xu){Xb — Xu)] , 



On the other hand, 

(/ + 2) 2a;Ta;H + {xh - xw + a;„)^ 



9f 



9l 



9l 



9l 



J{(/ + 2)(2 
+/ I^SxtXw + {xt + 3x„) 



Xh - 2xw + a;* + a;„)^ - 2 



2 2 
95 =96 



^XT + (/ + 2)(a;t + 3a;„) (2a;if - 2a;M^ +xt+ Xu) 
8a4 + ^{2xw + xt + 'iXu)xT + {xt + 3x„)^ | , 
(/ + 2) 2xTXH + {xh -xw+xt)^ , 

i{(/ + 2)(2a:H - 2xw + Xt + Xu)^ - 2 4xt + {f + 2){xu + 3xt) {2xh - 2xw + Xt + Xu) 
+/[8a;Ta;w. + (a;„ + 3xt)^ +2 8a;| + 4:{2xw + Xu+ 3xt)xT + (a;„ + 3a;t)^] |, 
i|(/ + 2)(2a;H - 2xvi/ + + a:„)^ + 2 /(a;T - a;t - a;„) - 2(a;t + a:„) {2xh - 2xw + Xt + 
+8xTXw{f + 2) + 8a;| - 2f{xt + Xu)xt + (/ + 2) (a;* + a;„)2|. 



(120) 



(121) 

(122) 
(123) 
(124) 
(125) 
(126) 
(127) 



(128) 

(129) 
(130) 

(131) 

(132) 
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C. The process 77 — > A<j)a 



Bosonic loops do not contribute to this process. The fermionic contributions are given through diagrams of the 
type (i) in Figs. [T] and [31 The corresponding amphtudes can be written as fohows: 



Ai — Fii + F13, 

A-2 = F23, 



(133) 
(134) 



where, the subscript j in Fij stands for the contribution of the particular set of diagrams. Using the same notation 
defined above, the partial amplitudes read 



Fii^F,i+Ff„ 

where F^ (Fii) represents the contribution due to Higgs boson A {Z boson). We can write these factors as 



f 
1 



sCo(l,2), 



327r3c^ 1-XZ+ iy/x^lz 



± L — Xl 

Gza.p^'^QaQ} J dxi J dx2 



X1X2 



Xf - X1X2 



(135) 

(136) 
(137) 



The expression given in Eq. (|137p corresponds to the one given in Ref. [67[ . 

^13 = Y.^fQ}{^AffQ4>Jf)[^)-^[^T{xA~ Xa)sD^{l,i.2) + 



2xf{l + XA- xa)s [Da{l,2, 3) + i?o(2, 1, 3) + Do{l, 3, 2)J + 

{Xa - Xa) [{Xt - Xa)Co% 3) + {xt " Xa)Co(1, 4) + (x„ - Xa)Co{l. 3) + {Xu " Xa)Co{2, 4)] }, (138) 



= i:^fQ){QAnQ,jf){^)s{s[x\^^^^ 

2xf{xu - xt)s [Do{1, 2, 3) + Do{2, 1, 3) + Do{l, 3, 2) 
I^ XaXA-xl \ r _ ^^)Coil, 3) + (x„ - xa)Co{2,4) 

\ Xt / l 



Xt 

Xa ~ Xf 



^Xa_ 

ixu-xt)[[^^±^)Co{1.2) + [ 



(xt - Xa)Co(2, 3) + [xt - xa)Co{IA) 

AXaXA - {xt + Xu)^ 



XT 



o(3,4)]} 



(139) 



V. NUMERICAL RESULTS AND DISCUSSION 



We now turn to discuss our results. We refer to our three benchmark scenarios already discussed, namely, (1) 
Scenario I (the decoupling limit), (2) Scenario II ( SM-likc), and (3) Scenario III (a more general case of 2HDM-III). 



A. Scenario I (the decoupHng limit) 

The main purposes of this subsection is to show explicitly how the decoupling of the charged Higgs effects operate 
according the criteria established in Ref. [HO]- In this scenario, one assumes the existence of a light Higgs boson h, 
with mass of order of the Fermi scale v, whereas the rest of Higgs bosons are assumed very heavy, i.e., irLfj± ~ niA ^ 
nfiH » V. The heavy Higgs boson effects decouple through two essentially different mechanisms [50[, namely, by 
assuming /i^2 or /i^2 but taking large tan/3 or cot /3, depending on the configuration chosen for the Ag 

and Ay parameters. Accordingly, we will consider the following two cases: 
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Case A m/i = 120 GeV, a = /3 — 7r/2, Ag = A7 = 0, fif^ » w^, and >> mh- 
Case B TO/i = 120 GeF, a = /? - 7r/2, Ag -A7 = 0.1, ^1^2 ^ "^^^ and >> m/j. 



Ca5e .4 

SM - 

(Jiz = 400 GcV, tani5= 1 - 
(JH = 400 GeF, taii/3=5 - 
/Ji2 = 500 GeV, taE/3=l - 
/<,2 = 500 GeV, tan/3= 5 - 



1 . \ . 

\ ■, \ \ 
\ \ \ 

\ '■, ^ 
\ '. \ 
\ '■, 



Case B 

SM - 

JUA = 900 GeV, tan /3 ^ 50 - 
mj = 1000 GeF, taa/3=50- 
mj = 1000 GeV, tan/3 = 80 - 
m,4 = 1100GeV; tan,3 = 80- 



500 600 700 800 900 500 600 700 800 900 

mffi {GeV\ 

FIG. 6: Behavior of cross section for the process 77 — >■ /i/i as a function of the charged Higgs mass in the Scenario I (decoupling 
hmit). The diverse values of the parameters as well as the case in consideration (A or B) are shown in the figure. In all cases 
the center-of-mass energy -^s = 500 GeV was used. 

In the Figure m we present the cross section of the process cr(77 — hh) as a function of the charged Higgs mass rn/fi . 
One can get the SM result (represented in these figures by the horizontal thick line) when m^fi >> m/t, showing the 
decoupling nature of this contribution. On the other hand, we can see from Fig. [7] that in the limit tan/3 << 1 of 
the case A the SM result is obtained. The decoupling nature of the contributions for the case B is clearly appreciated 
when tan/3 >> 1. 

10^ [ — I 1 1 1 1 1 1 1 ] "1 I 1 1 1 1 1 1 1 



P12 = 400 GeV, 41 r 300 GeV 

SM 

mfft^SOO GeV 

mH' = 600 GeV 

m.H± = 700 GeV... 



771^* = 500 GeV 

m^i = 600 GeV 

m». = 700 GcV- - -- 
.1 = 850 GeV = 500 GeV 

Case B 



tan/3 



FIG. 7: Behavior of the cross section for the process 77 hh as & function of tan/3 in the Scenario I (decoupling limit). 
The diverse values of the parameters as well as the case in consideration (A or B) are shown in the figure. In all cases the 
center-of-mass energy -^s = 500 GeV was used. 
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B. Scenario II ( SM-like) 

As already commented, in this scenario the couphngs hVV {V — W,Z), hhh, hhhh are nearly indistinguishable 
from the corresponding ones of the SM, but the hff couplings can deviate significantly from their SM counterparts 
hsnff- In this context, we will use the values nifi = 120 GeV, niA = 110 GeV, niH = ra^i — mA + m^, and 
Mi2 = 130 GeV in the two cases (Ae = —Ay — 0.1) and (Ae = — A7 ~ 1). The values Xuu,dd ~ 1, —1, which arise from 
a selection of a specific texture of the Yukawa matrices, will be used. In addition, we will assume that a — f5 ± tt/2. 




300 350 400 450 300 350 400 450 

-ys [GeV\ 

FIG. 8: Behavior of cross section for the process 77 /i/i as a function of center-of-mass energy ^/s in the Scenario II (SM-like). 
The curves correspond to SM (solid line), 2HDM type I-II (dashed line), and 2HDM-III (dashed-dotted line), for tan/3 — 5 
(left), 15 (right) in the case Ae = — A7 = 0.1 (up panels), 1 (down panels). 

In Fig. [SJ the cross section for the process 77 —> /i/i as a function of center-of-mass energy ^/s is shown. It can be 
appreciated that for tan /3 = 5 both the SM and the 2IIDM's predictions essentially coincide. However, one can see 
that for tan/3 = 15 the cross section predicted by the 2IIDMs could be two order of magnitude larger than the SM 
result. The results for x = 1, —1 are very similar. In this case, the cross section of the mode 17(77 ~^ f^f^) ^ "^0 fb. 

On the other hand, in Fig. ^ the cross sections for the 77 — >■ AA reaction as a function of the center-of-mass 
energy -^/s is displayed and the predictions of the 2HDM-III compared with those generated by the 2HDM-II. In this 
process we also consider the contribution of the parameters of the Higgs potential (Ae = — A7 = 0.1) (up panels) and 
(Ae = — A7 = 1) (down panels). We can observe that the main impact for the cross section comes from (Ag = — A7 = 1) 
and tan/3 large. The cross section of this mode could be enhanced by two orders of magnitude compared with the 
case Ag = — A7 = 0.1 or the usual case Xq = -A7 = (2IIDM-II). One can get (7(77 — >■ AA) ~ 5 x 10^ fb for 
tan/3 = 15, with y/s — 450 GeV, taking x = 1, — 1 and Ag = — A7 = 1. The results for the cross section of the 
process 77 — > HH are shown in Fig. [10] for the same parameters of the previous process. Likewise, the cross section 
predicted by the 2HDM-III is two orders of magnitude larger than the one predictied by the 2IIDM-II. We can obtain 
cr(77 HH) ^ 1 X 10^ fb for tan/3 = 15, x = Ij —1; with ^/s = 500 GeV and Ag = — A7 = 1. As far as the cross 
section for the (7(77 — > hH) process is concerned, its behavior as a function of the center-mass energy ^/s is shown 
in Fig. [TT] From this figure, a considerable enhancement of the cross section compared with the 2HDM-II prediction 
can be observed for Ag = 1 = — A7. In this case, (7(77 — )■ hH) ~ 5 x 10^ fb for tan/3 = 15, around ^/s = 470 GeV and 
X = 1, — 1 . It is worth commenting that this is the first time that this process is studied. Another process which has 
not been studied in the literature is (7(77 — J> hA). The corresponding cross sections as a function of the center-mass 
energy ^/s is shown in Fig. [121 It can be appreciated from this figure the importance of the cross section for the case 
Ag = 1 = — A7 and tan/3 = 15, which is quite large as compared with the prediction of 2IIDM-II where the Ag and A7 
parameters are absent. It can be appreciated that the cross section could be of the order of 7 fb for y/s — 350 GeV. 
On the other hand, when Ag = — A7 << 1 the cross section is very insignificant to be considered as relevant signals 
of neutral Higgs bosons. The last numerical results of this scenario is the cross section for the process 77 — >■ HA, 
which is shown as a function of the center-mass energy in Fig. 1131 It can appreciated an important value for the cross 
section of about 2 x lO'^ fb for the case Ag = — A7 = 1, x = ^ li tan /3 = 15, and ^/s around 350 GeV. 
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FIG. 9: Behavior of the cross section for the process 77 — >■ AA as a function of the center-of-mass energy ^/s in the Scenario 
II (SM-hke). The hnes correspond to 2HDM-II (sohd hne), 2HDM-III with Xnu,dd = 1 (dashed hne) and 2HDM-III when 
Xuu,dd = —1 (dashed-dotted hne), for tan/3 = 5 (left), 15 (right) in the cases Ae = — A7 = 0.1 (up panels) and Ae = — A7 = 1 
(down panels). 
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FIG. 10: Behavior of the cross section for the process 77 — >■ HH as a function of the center-of-mass energy ^/s in the Scenario 
II (SM-like). The description of the plots is the same as in Fig. [9] 

C. Scenario III (a more general case of 2HDM-III) 

As already commented, this scenario is much more general than the scenario II, because arbitrary couplings of 
neutral Higgs bosons to SM particles are assumed. Also, the contributions of the Higgs potential Ag and A7 parameters, 
as well as the contributions of the Yukawa texture in the couplings (t>ff, are included. As commented at the end of 
Sec. nil the degenerate and the nondegenerate cases, as well as the case with a light CP-odd scalar will be considered. 
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FIG. 11: Behavior of the cross section for the process 77 — >■ hH as a function of the center-of-mass energy in the Scenario 
II (SM-Hke). The description of the plots is the same as in Fig. (9] 




300 400 500 600 700 800 900 '" 300 400 500 600 700 800 900 



FIG. 12: Behavior of the cross section for the process 77 — > hA as a function of the center-of-mass energy ^/s in the Scenario 
II (SM-like). The lines correspond to 2HDM-II (solid line), 2HDM-III with Xv.u,dd = 1 (dashed line) and 2HDM-III when 
Xuu,dd = —1 (dashed-dotted line), for tan/3 — 5 (left), 15 (right) in the case Ag ~ —Xj = 1. 



1. The nondegenerate case 

We first discuss the nondegenerate case, defined by the values mjy± = 400 GeV, tua = 350 GeV, Mh = 520 GeV, 
fii2 = 120 GeV, and rrih = 120 GeV. The set of values Ay = —Xq = —1 and A7 — — Ag — —0.1 are considered. In 
addition, it is assumed that a = /3 and a — /3 ± 11/2. 

In the down panels of Fig. [Ml we show the cross section for the process 77 /i/i as a function of the center-of-mass 
energy ^/s. An important difference can be appreciated between two specific values of the Yukawa matrices with 
textures x = 1 and % = — 1. For tan {3 — 5 and x — the cross section could be up to one order of magnitude larger 
than for the case x = 1, when ^/s < 350 GeV. In this region, we can get (7(77 — hh) ~ 1 x 10^ fb for -^s around 
370 GeV. The cross section predicted by the 2HDM-II is two (one) order of magnitude lower than the 2HDM-III 
prediction with x = ~l(l)- It can be seen that for -^i > 500 GeV and x = 1 ^nd x — the cross section is of 
the same order of magnitude. However, in the 2HDM-III the cross section could be larger than the result obtained 
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FIG. 13: Behavior of the cross section for the process 77 — >■ HA as a function of the center-of-mass energy in the Scenario 
II (SM-like). The description of the plots is the same as in Figure (9] 

in the context of the 2HDM-II, due to the following choice Xq — —Ay = 1 of the parameters. In the same figure, 
it can be appreciated a spectacular enhancement of the cross section of cr(77 — > hh) ~ 3 x 10^ fb for tan/3 — 30, 
X = — 1, and y/s = 350 GeV. The cross section predicted by the 2HDM-II is one order of magnitude lower than the 
one predicted by the 2HDM-III. On the other hand, it can be appreciated from the up panels of the same figure that 
in Ag = — A7 = 0.1 case, the corresponding cross sections are suppressed by about one order of magnitude with respect 
to those obtained in the Ae = — A7 = 1 case. 

On the other hand, the cross section for the 77 — > HH reaction is shown in Fig. [15] as a function of the center-of- 
mass energy ^/s. The signal for this process could be relevant in the TeVs region, y/s > 1 TeV. Therefore, this mode 
could be far away of the reach of early linear colliders. It can be appreciated from this figure that the predictions of 
the THDM-III approximates to that of the THDM-II in the case Xg = -A7 = 0.1. 

As far as the 77 — >■ AA process is concerned, the corresponding cross section as a function of the center-of-mass 
energy y/s is shown in Fig. [TH This mode is important in the 2HDM-III for center-of-mass energies above 700 GeV 
and for large values of tan/3. From this figure, it can be appreciated that in the case Ag = — Ay = 1 (up set of graphics 
inf Fig. [TH), cr(77 ^ AA) - 1 x 10^ fb for = 800 GeV, tan/3 = 30, and a = /3 or a ^ ± tt/2. The 2HDM-III 
predictions are two orders of magnitude larger than the ones of 2HDM-II. In the scenario with tan/3 = 5, the cross 
section predicted by 2IIDM-III is of order of 10 fb for ^/s = 800, x = 1, and a — (3 , whereas the 2HDM-II prediction 
is about one order of magnitude lower. However, the situation changes drastically when a = /3 ± 7r/2, as in this case 
the 2HDM-II contribution dominates. On the other hand, it can be appreciated from these figures that in the case 
Ag = — Ay = 0.1 (down set of graphics in Fig. [T5| . the corresponding cross sections are of the same order of magnitude 
that those predicted by the THDM-II. 

We now turn to discuss the process 77 — > hH. The corresponding cross section as a function of the center-of-mass 
energy -^s is shown in Fig. I17[ in which the up set of graphics corresponds to the case Ag = —Ay = 1, whereas the 
down set was obtained using the values Ag = — Ay = 0.1. It is found that this process is sensitive to tan/3 and the 
mixing angle a. The cross section can reach a value of 20 fb for tan/3 = 5 and a = /3. For large tan/3 values, the 
cross section is enhanced by several orders of magnitude. In fact, cr(77 — hH) ~ 1(5) x lO** fb for y/s = 800 GeV, 
tan/3 = 30, a = /3 (a = /3 ± 7r/2), and x — ±1- From these figures, it can be appreciated that the prediction of the 
THDM-III for the Ag = —Ay = 0.1 case is about one order of magnitude lower than that for the Ag = —Ay = 1 case, 
and clearly tends to the THDM-II prediction. 

The cross section for the 77 — > hA as a function of the center-of-mass energy ^/s is shown in Fig. [181 in which the up 
set of graphics corresponds to the Ag = —Ay = 1 case, whereas the down set of figures arises from the Ag = —Ay — 0.1 
case. This cross section is quite sensitive to the mixing angle a. A relevant value for the cross section arises when 
tan /3 is large and a — /3. In fact, cr(77 — > hA) ~ 1 x 10* fb for y/s = 500 GeV, tan /3 = 30, and a = /3. It is important 
to notice that in the case Ag = —Ay = 0.1, only results for a = (3 are presented, as for a = /3 ± 7r/2 the cross sections 
are essentially independent on the Ag and Ay parameters, being therefore almost identical to those of the down panel 
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FIG. 14: Behavior of cross section for the process 77 — > /i/i as a function of the center-of-mass energy ^/s in the Scenario 
III (the nondegenerate case). The lines correspond to 2HDM-1I (sohd hne), 2HDM-III with Xuu.dd ~ 1 (dashed hne) and 
2IIDM-1I1 when Xv,u,dd = —1 (dashed-dotted line), for tan/3 = 5 (left), 30 (right) in the cases Ae = — A7 = 0.1 (up panels) and 
Ae — — A7 = 1 (down panels). We choice the mixing angle a = /?. 

of the up set of graphics. It can be appreciated from these figures that the cross sections for the Ag = —Ay = 0.1 case 
are about one order of magnitude lower than those for the Ag — —Ay = 1 case. 

The cross section for the 77 — >■ HA process is shown in Fig. [TO] as a function of the center-mass energy, in which 
the up set of graphics corresponds to the Ag = — A7 = 1 case, whereas the down set of figures arises from assuming 
Ag — — A7 — 0.1. Relevant cross sections are predicted by the 2HDM-III of order 10^ fb for tan/3 = 30, x = ^^^d 
energies around 900 GeV. It is important to notice that in the case Ag — — A7 = 0.1, only results for a = /3 + 7r/2 are 
presented, as for a — (3 the cross sections are essentially independent on the Ag and A7 parameters, being therefore 
almost identical to those of the down panel of the up set of graphics. It can be appreciated from this figure that for 
tan/3 = 5 the cross section in the Ag = —Ay = 0.1 case is one order of magnitude larger than in the Ag = —Ay = 0.1 
case, in contrast with the behavior observed in all the other processes. 

2. The degenerate case 

In this paragraph, we present results for the degenerate case, which is defined in Sec. II. Only the case Ag = —Ay = 1 
will be considered, as the Ag = —Ay = 0.1 case leads to cross sections suppressed by about one order of magnitude 
with respect to the former one. Although in general terms the cross sections for tan/3 = 30 tend to be about two 
orders of magnitude larger than those obtained using tan/3 = 5, we have preferred to make predictions using only the 
latter value because in this case the predictions of the two versions of the model (THDM-III and THDM-II) can be 
clearly distinguished. 

In Fig. [201 the behavior of the cross sections for the processes 77 — ?> hh and 77 — > HH are shown as a function 
of the center-mass-energy, with a — jS in the former process and a = /3 ± 7r/2 in the latter one. Besides to optimize 
the cross sections, these choice of values maximize the differences between both models. It can be appreciated from 
this figures that the THDM-III predicts cross sections as large as 10^ fb and 10 fb for the hh and HH channels, 
respectively, which are two and one orders of magnitude larger than those predicted by the THDM-II. 

The cross sections for the processes 77 — ?► AA and 77 — > hH are shown in Fig. [21] as a functions of the center- 
mass-energy, ior a — f3 in both cases. From this figure, it can be appreciated that the THDM-III prediction for the 
cross sections of both processes range from about 1 fb to 10 fb in the energies range shown. In contrast, the THDM-II 
predict cross sections quite suppressed (10^^ fb for 77 — > AA and 10^^ fb for 77 — >• hH), which varies slightly in all 
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FIG. 15: Behavior of cross section for the process 77 — >■ HH as a function of the center-of-mass energy ^fs in the Scenario III, 
with a nondegenerate masses spectrum. The hues correspond to 2HDM-II (sohd hne), 2HDM-III with Xuu.drf = 1 (dashed line) 
and 2HDM-III when Xtm.dd = ~1 (dashed-dotted hne), for tan/3 = 5 (left), 30 (right) in the cases As = — A7 = 0.1 (up panels) 
and Ae = — A7 = 1 (down panels). We choice the mixing angle a = /3 ± 7r/2. 



the energies range considered. 

The cross sections for the processes 77 — )■ hA and 77 — HA are shown in Fig. [52] as a function of the center- 
mass-energy. The relations used between the a and /3 angles are shown in the figure. In this case, the predictions of 
the THDM-III ranges from 1 jh to 10 /& for the 77 — hA process, whereas the prediction for the 77 — ^ HA reactions 
is one order of magnitude lower. The predictions of the THDM-II are quite suppressed, as both cross sections are of 
order of 10"^ /6 or lower. 



3. A light CP-odd scalar 



In this paragraph, we discuss the very interesting case of a light CP-odd scalar A, which is allowed by the current 
constraints on the model. We consider the three possible processes, namely, 77 — > AA^ 77 — > HA, and 77 — > hA. In 
Fig. [23l the behavior of the cross sections for these processes as a functions of the center-mass-energy is shown, in a 
scenario with tjia = 50 GeV, mi, = 120 GeV, mH± = 350 GeV, m// = 400 GeV, fii2 = 70 GeV, and tan/3 = 5. The 
value a = (3±tt/2 for the 77 — > AA and 77 HA processes is assumed, whereas in the case of the 77 — >■ hA reaction 
it is assumed that a — /3. In all these processes, it is assumed that Ag = —Ay = 1 and that Xuu — Xdd = {1, — !}• 
Our notation and conventions are shown in the first graph of Fig. 1231 From the first graph of this figure, it 
can be appreciated three resonant effects for the 77 — AA process, centered at energies ^/s = 120 GeV — mh, 
\/s = AGO GeV = mn, and ^/s — 350 GeF = 2mt. The resonant effects due to nih and mn are spectacular, as the 
cross section can reach values of up to 10^ fb and 10^ fb, respectively. The resonant effect at 2mt is less significative, as 
it occurs through a l~loop fluctuation. Apart from these resonant effects, the values of the cross section are within the 
range of variation encountered in other scenarios analyzed previously, as the cross section predicted by the THDM-III 
ranges approximately from 10^^ fb to 1 fb in both the Xuu = Xdd = — 1 and Xuu = Xdd = 1 scenarios, whereas in the 
THDM-II the corresponding cross section is about one order of magnitude larger. As far as the 77 — > HA process 
is concerned, it can be seen from this figure that the cross section predicted by the THDM-III ranges from 10~^ fb 
to 1/6 for the scenario with Xuu = Xdd = — 1, whereas for Xuu = Xdd = 1 the cross section ranges from 10^"^ /6 
to 10~^ fb. The corresponding cross section predicted by the THDM-II ranges from 10~^ fb to 10~^ fb in the same 
range of variation of the center-mass-energy. Finally, it can be appreciated from the third graph of Fig. [23] that the 
cross section for the 77 — ?► hA reaction ranges from 10 fb to 10^ fb in the THDM-III in the scenario %„„ = Xdd = — 1, 
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FIG. 16: Behavior of cross section for the process 77 — > A A as a function of the center-of-mass energy -^s in the Scenario 
III, with a nondegenerate masses spectrum. The lines correspond to 2HDM-II (solid line), 2HDM-III with Xuu,dd ~ 1 (dashed 
line) and 2HDM-III when Xnu,dd = —1 (dashed-dotted line), for tan/3 = 5 (left), 30 (right), taking a — 13 (up panels) and 
a = P ± ■k/2 (down panels). The up(down) set of graphics corresponds to the case Ae = — A7 — 1 (Ae — — A7 = 0.1). 



and from 1/6 to 10 /& in the scenario Xuu = Xdd = ^1- The corresponding cross section in the THDM-II is quite 
suppressed, as it ranges from 10~^ fb to 10~^ fb in the same domain of energies. The cross sections for all these 
processes can be enhanced by at least in one order of magnitude if tan/3 = 30 is used instead of tan/3 = 5. 
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FIG. 17: The figure shows the behavior of cross section for the process 77 — >■ hH as a function of the center-of-mass energy 
in the Scenario III for the nondegenerate case. The description of the plots is the same as in Figure [T^ 



VI. SUMMARY 

In this paper, a comprehensive study of the one-loop 77 — >■ 4>i(l)j {(j^i = h, H, A) processes in the context of a general 
version of the two Higgs doublet model (2HDM-III) was presented. A nonlinear i?^-gauge, which allows us to define 
the W gauge boson propagator in a covariant way under the electromagnetic gauge group, was used. This gauge, 
which reduces significantly the number of Feynman diagrams to be considered in comparison with those that must 
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FIG. 18: The figure shows the behavior of cross section for the process 77 — >■ hA as a function of the center-of-mass energy 
y^s in the Scenario III for the nondegenerate case. The description of the plots is the same as in Figure [161 In the down set of 
diagrams, only the case a = /3 is considered. 

be calculated in conventional linear gauges, makes the issues of electromagnetic gauge invariance and cancelation of 
divergencies simpler. Explicit analytical expressions for all the possible modes, namely hh, HH, AA, hH, hA, and 
HA, were presented. The corresponding amplitudes are completely general in the sense that they can be used for 
any version of the 2HDM. The version of the 2HDM that is considered in this work, which we called simply 2HDM- 
III, comprise the implementation of a flavor symmetry in the Yukawa sector, namely a four-zero Yukawa Texture, 
which allows us to suppress FCNC effects without necessity of introducing a discrete symmetry. Due to this, a Higgs 
potential, more general than the one considered in the 2HDM-II version, can be introduced. This Higgs potential 
includes two dimensionless parameters, Xq and A7, to which the cross sections for the 77 — >■ (f>i(j)j processes are quite 
sensitive. The 77 <j)i4>j mechanisms for Higgs pair production were analyzed in three scenarios of the 2HDM-III. 
In the scenario I (the decoupling case), the 77 — ^ hh reaction was studied in the two possible cases in which the 
decoupling operates, namely, when /i^2 ^ or 1A2 ^ but assuming that tan/3 or cot/3 are large, depending 
on the configuration chosen for the Ag and A7 parameters. In both cases, it is found that, in the heavy mass limit 
of the charged Higgs, the cross section for this reaction approach to the well known SM result. In the scenario II 
(SM-like), it is assumed that the hVV {V = W,Z), hhh, and hhhh couplings are nearly indistinguishable from the 
corresponding ones of the SM, but the hff couplings can deviate significantly from their SM counterparts hsAiff - It 
was found that some combinations of Yukawa textures with the Ag and A7 parameters, together with large tan/3 lead 
to experimentally interesting cross sections. In the scenario HI (a more general case of 2HDM-III), besides considering 
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FIG. 19: Behavior of cross section for the process 77 — ^ HA as a function of the center-of-mass energy ^/s in the Scenario III 
for the nondegenerate case. The description of the plots is the same as in Figure [T5] In the down set of diagrams, only the 
case a = /3±7r/2is considered. 



the contributions of the Higgs potential through the Xq and A7 parameters, as well as the contributions of the Yukawa 
texture in the couplings 0//, arbitrary couplings of neutral Higgs bosons to SM particles were assumed. In this 
scenario, the implications of a light CP-odd scalar was studied. 

In general terms, we can conclude that the parameters of the Higgs potential Ag and A7 considerably enhance the 
cross sections for the 77 — >■ processes. In almost all cases, the results of 2HDM-III are two orders of magnitude 
larger than those obtained from the 2HDM-II. A considerable enhancement for the cross sections of the processes 
77 — >■ (j)i4>j was observed in the regime of large tan /?. 



Appendix A: Feynman rules in the nonlinear gauge 

We present the couplings for the most polular types of 2HDM (I, II, y HI). The vertex can be written as 

where G^-ff are dimensionless functions expressed in the table (jl]). 



1. Yang-Mills Couplings 

This sector is strongly affected by the gauge-fixing procedure that we used. This couplings can be found in (4^ . 
Here, we only present the Feynman rules associated with the vertices used in this work, namely, Ajj{ki)W^ {k2)W~ {k^) 
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FIG. 20: Behavior of cross section for the processes 77 — >■ hh (left) and 77 — > HH (right) as a function of the center-of-mass 
energy -^s in the Scenario III for the degenerate case. 
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FIG. 21: Behavior of cross section for the processes 77 — >■ AA (left) and 77 — > hH (right) as a function of the center-of-mass 
energy ^/s in the Scenario III for the degenerate case. 



and AaApW^Wp . The corresponding vertex functions are given by 



ieTxp'q{ki,k2,k:i) = -ie ( (fca - A;2)i,5Ap + (fci ~ ^3 - ^k2)\gpr, + {k2 - ki + ^k3)pgx,j ) , (A2) 



ie^'^apxp = -ie'' { ^g^^gxp - (1 - ^){ga\g0p + gapgpx) 



(A3) 
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FIG. 22: Behavior of cross section for the processes 77 — >■ hA (left) and 77 — >■ HA (right) as a function of the center-of-mass 
energy ^/s in the Scenario III for the degenerate case. 
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FIG. 23: Behavior of cross section for the processes 77 — > AA, 77 — > HA, and 77 — > HA as a function of the center-of-mass 
energy in the scenario of a hght CP-odd scalar. 
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Qhuu 

GHlliQudd 



Type-I 



Type-II 



Type-III 



+ 



-fiXdd 



+ 



V2c 

C-I3XV.V. 

a~l3Xdd 
V2cg 



V2s 



Qau,Qa 



-tp + 



QAuu 



^/2s 



TABLE I: Dimensionless function that define the Yukawa couphngs [13, [28l . [30 



Scalar and Kinetic sector 



Now we present the couplings that contain scalar particles, in particular the vertex functions necessary for our 
calculation. The vertex functions between scalar Higgs can be written as following: 



9<piH±HT 

9h±g1.a 

9<l>aAA 



4 



Scj>a<Pb4'i: I 



16 



■Qh± 



± 



'igmw , 



AA 



(A4) 
(A5) 
(A6) 
(A7) 
(A8) 



The explicit form of this vertex function can be found in [50|. 

The couplings that involve a pair of of VK-Goldstone boson was modified by the gauge-fixing procedure. This 
couplings can be expressed as: 



GipiWw, 



16 ^G±G^0i0j 



(A9) 
(AlO) 



The explicit form of this couplings have to be take of the reference [44 1. 

On the other hand, the couplings involving gauge and scalar fields are present in the tableHIl Here, the dimensionless 
functions Q can be written as: 



QhWW = —Gw±H^H = —GzAH =S(^_q), 
GhWW = Qw±H^h = GzAh — C(^p~a}- 



(All) 
(A12) 



Appendix B: The scalar functions 



In this appendix we present the Co(a) and Do{a) scalar functions appearing in the amplitudes Ai^ and A25 associated 
with the processes 77 — > A A and 77 — > (t>a4'b- 
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TABLE II: Structure of the couplings involving gauge and scalar fields. 



Coupling 


Vertex Function 


Coupling 


Vertex Function 


ffAf,H-(p_)H+{p+) 




9a^a^h- H+ 




^^/i<3w(P-)G5,(p+) 


je(p_ — p+),, 


i^A^A^G^G^ 








9(hh,HH,AA)W^W^ 




S'0o(Pa)VK±HT(p„±) 




9A^Wi^H^(l>a 




^A(p^)W±iifT(p^±) 




9Ai_,W^H^A 








gZA(pA)<t>a{pa) 


2^(Pa - PA)^,QzA^a 



Co(l) = Co(s,m-,m|,TO^,m^,m^±), (Bl) 

Co(2) = Co(s,m-,m|,m^±,m^±,m^), (B2) 

Co(3) = Co(0,0,s,m^±,m^±,m^±), (B3) 

Co(4) =Co(0,0,s,m2^,m2^,m2^), (B4) 

Co(5) = Co(0, mj,u, m^±,TO^±,TO^), (B5) 

Co(6) = Co(0, mj,u, m^,m^,m^±), (B6) 

Co (7) = Co(0, m-,i,m^±,m^±,m^), (B7) 

Co(8) = Co(0,mi,t,rn^,m^,m^±), (B8) 

Co(9) =Co(0,^,m2,m^±,m^±,m2^,), (B9) 

Co{lO) = Co{0,t,mlm^^,ml,,mj,±), (BIO) 

Co(ll) = Co(0, u,m|,m^±,m^±,m^), (Bll) 

Co(12) = Co(0,u,m2,m2^,m^,m^±), (B12) 

I?o(l) = £)o(0, 0,m^,m|,s, u,m^±,m^±,m^±,m^), (B13) 

i:>o(2) = L'olO, 0,TO^,rn|,s, u,TO^,m^,TO^,m^±), (B14) 

i:>o(3) = i'olO, 0, TO^, to|, s, TO^± , TO^± , TO^± , rn^), (B15) 

£'0(4) = L'olO, 0,m^,m|,s,t,m^,m^,m^,m^±), (B16) 

I^oCS) = 1)0(0, mj^,0,m^,u,i,m^±,m^±,m^,m^), (B17) 

£'0(6) = 1)0(0, m,^,0,m^,u,i,m^,m^,m^±,m^±), (B18) 

where m,; = mj = niA for 77 — > ylA, and = ma and mj ~ mi, for 77 — >■ (jia't'b- 
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